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Abstract 

The coupling of a dilaton to the 5'J7(2)-Yang-Mills field leads to interesting non-perturbative 
static spherically symmetric solutions which are studied by mixed analitical and numerical meth- 
ods. In the abelian sector of the theory there are finite-energy magnetic and electric monopole 
solutions which saturate the Bogomol'nyi bound. In the nonabelian sector there exist a count- 
able family of globally regular solutions which are purely magnetic but have zero Yang-Mills 
magnetic charge. Their discrete spectrum of energies is bounded from above by the energy of 
the abelian magnetic monopole with unit magnetic charge. The stability analysis demonstrates 
that the solutions are saddle points of the energy functional with increasing number of unstable 
modes. The existence and instability of these solutions are "explained" by the Morse-theory 
argument recently proposed by Sudarsky and Wald. 
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1 Introduction 

As it is well known, the Yang-Mills (YM) equations arc scale invariant which excludes glob- 
ally regular (i.e. non-singular with finite energy) static solutions [1,2]. The usual method for 
circumventing this nonexistence result is to introduce a Higgs field. The coupling of the Higgs 
field has two effects. First, it breaks scale invariance. Second, a nonabelian gauge group G gets 
spontanously broken to a subgroup H. If the homotopy group TTD-iiG/ H) is nontrivial (where 
D is the number of space dimensions), then the coupled YM-Higgs theory has topologically 
stable solutions. A prominent example is the t'Hooft-Polyakov monopole [3] in the SU {2)-YM. 
theory with a triplet Higgs field. 

A spontanously broken gauge theory may admit another class of globally regular solutions 
if itd(G/H) is nontrivial. This homotopy group is isomorphic to the group of loops in the 
configuration space (i.e. space of static, finite energy configurations). Nontriviality oi TToiG/ H) 
means that there are noncontractible loops passing through the vacuum. The argument, due to 
Taubes [4] and Manton [5] , of how such noncontractible loops lead to a nontrivial solution runs 
as follows. Consider all loops starting and ending at the vacuum in a fixed homotopy class. On 
each loop there is a configuration of maximal energy and the infimum of these energies gives 
a saddle point of the energy functional (and therefore the static solution). Due to the non- 
compactness and infinitc-dimcnsinality of the configuration space, this argument is obviously 
not rigorous, and to actually prove that the mini-max procedure converges is a difficult technical 
problem. Static solutions corresponding to saddle points of the energy functional were called 
sphalerons to emphasize that, in contrast to solitons, they are unstable. The existence of a 
sphaleron was first shown by Taubes [4] in the SU{2)-YM theory with a triplet Higgs field and 
by Manton [5] in the SU{2)-YM theory with a complex doublet Higgs field. 

Although sphalerons were originally discovered in spontanously broken gauge theories, it 
should be stressed that the Higgs mechanism is by no means necessary for the existence of 
a sphaleron. Actually, this is already clear in the SU {2)-YM theory with a complex doublet 
Higgs, where the gauge group SU (2) is completely broken and the homotopy group relevant for 
constructing a sphaleron is ■K3{SU (2)) ~ Z. Thus, in this case the role of the Higgs field is just 
to break the scale invariance while the gauge group itself has a nontrivial third homotopy group. 
This suggests that a sphaleron may exist in the SU{2)-YM theory coupled to other fields (of 
attractive force), provided that the coupling: i) breaks scale invariance, and ii) does not alter 
the topology of the configuration space of pure 5f7(2)-YM theory. 

In this paper I consider a simple example of a coupling which satisfies these two requirements, 
namely the coupling of a dilaton. The dilaton, (p, is a real (massless) scalar field which couples 
to other matter fields (with lagrangian Lm) through the term e^^"'^L„i, where a is the dilaton 
coupling constant. I will show that static spherically symmetric S'C/(2)-YM-dilaton equations 
have globally regular solutions with the following properties: 

a) there exist a countable family of solutions Xn (n e N), 

b) the energy increases with n and is bounded from above, 

c) the solution X„ has exactly n unstable modes. 

This family of solutions is in striking analogy to the Bartnik-Mckinnon (BM) solutions [6] of the 
Einstein- S*?/ (2)-YM equations, which have the same properties a)-c). In both cases the solution 
Xi may be interpreted as a sphaleron (for BM solutions this was first observed by Mazur [7]; 
see also [8]). 
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A natural question arises: why do two theories with completely different dynamics have 
qualitatively the same spectrum of solutions? The answer was recently proposed by Sudarsky 
and Wald (SW) [9]. They presented a heuristic but convincing argument which accounts for 
the properties a)-c) (except for the boundcdncss of energy) in the case of BM solutions. This 
argument is formulated in the spirit of Morse theory for Hamiltonian systems and exploits 
the existence of topologically nonequivalent multiple vacua in the SU (2)-YM theory (which is 
related to the fact that 'it^{SU{2)) ~ Z. A detailed description of the SW argument will be given 
in Section 7. Here, let me only note that in the case of the solution Xi, the SW argument is, in 
essence, equivalent to the mini-max procedure for noncontractible loops. However, in contrast 
to the mini-max construction, the SW argument can be naturally extended (admittedly, under 
additional assumptions) to account also for the existence of solutions with n > 1. Although 
the SW argument was originally formulated in the context of Einstein- YM theory, it is essentially 
insensitive to the concrete form of coupling, and applies almost without modifications to the 
YM-dilaton theory. In this sense, SW predicted the existence of solutions found in this paper. 
On the other hand, the results of this paper give further credence to the SW argument. 

The existence of the upper bound for the spectrum may be understood by considering 
the U{1) sector of the YM-dilaton theory. Surprisingly enough, there are finite energy abelian 
solutions which describe magnetic and electric point monopoles. The finiteness of energy is due 
to the regulating effect of the dilaton which weakens the short distance singularity. Morever, 
these solutions saturate the Bogomol'nyi bound (in the U{\) sector), hence their energies are 
equal to their charge. It turns out from numerics that the limiting solution X^o (whose energy 
bounds the spectrum from above) corresponds to the abelian magnetic monopole with unit 
magnetic charge. 

The YM-dilaton theory and the Einstein- YM theory may be embedded in a single Einstein- 
YM-dilaton theory governed by the action 



(1) reduces to the YM-dilaton theory. When a = 0, the action (1) becomes the Einstein- YM 
theory (plus trivial kinetic term for the scalar field). Finally, the case a = 1 corresponds to 
the low-energy string theory. It was shown by the author elsewhere [10] that the theory defined 
by the action (1) has static spherically symmetric (globally regular and black-hole) solutions 
with properties a)-c), for all values of a. This paper specializes to the limiting case a ^ oo. It 
seems instructive to consider this case separately, because it involves the essential features of 
the general case, but has an advantage of being considerably simpler, which allows to obtain 
some analitical estimates on the parameters of solutions. Also, the non-perturbative effect of 
the dilaton can be clearly seen in this model. 

The paper is organized as follows. In the next Section the field equations are derived and 
some scaling properties are discussed. In Section 3 the explicit abelian solutions are described. 
In Section 4 the a priori behaviour of globally regular solutions is obtained. In Section 5 the 
numerical results are presented and some qualitative properties of solutions are discussed. The 
deep analogy between these solutions and the BM solutions is emphasized. Section 6 is de- 
voted to the stability analysis. Finally, in Section 7 the SW argument is summarized and some 
possibilities of proving rigorously the existence of numerical solutions are suggested. 




(1) 



This theory is characterized by a dimensionless parameter a = a/\/G. When a ^ co, the action 
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2 Field equations 

The dynamics of the SU{2)-YM field coupled to a dilaton is defined by the action 



S 



(2) 



where F = dA + eA A A is the YM curvature of the SU{2) connection A and (f) is the dilaton. 
Hereafter, for convenience I put the coupling constants o = e = 1, which is equivalent to 
choosing a/e as the unit of length and 1/ae as the unit of energy. 
The field equations derived from (2) are 



D(e-^<P *F) = , 



(3) 
(4) 



where D is the SU{2) covariant derivative. 

I wish to find static spherically symmetric solutions to these equations that are globally 
regular i.e. non-singular and with finite energy. 

The most general spherically symmetric SU{2) connection has the form [11] 



A = ar^dt + br^dr + {wti + dT2)d'd + (cot 'dr^ + wt2 — dri) sinMip , 



(5) 



where a, b, w and d are functions of {r,t) and r, {i = 1,2,3) are generators of su{2) Lie 
algebra. Using the residual gauge freedom the radial gauge 6 = can be imposed. When the 
connection is static, i.e. a, w and d depend only on r, one can also set d = by a constant gauge 
transformation. Hence, the general static spherically symmetric SU (2) connection is described 
by two functions: the electric potential a{r) and the magnetic potential w{r). Now, I assume 
further that a = 0. Actually, this is not a restriction because one can show, following the 
argument given in [12], that there are no globally regular solutions with nonzero electric field. 
The purely magnetic YM curvature is 



F = w'ridr Ad-d + w'T2dr A sin-ddtp — (1 — w'^)Tsd'd A siwddtp , 



(6) 



where prime denotes differentiation with respect to r. For F given by this ansatz and for 
(f) = (f)(^r) , the equations (3) and (4) reduce to the following system 



1 



w 



0, 



{r^(j)y + 2e 



-2(f> 



,2 {l-W-'f 



2r2 



These equations may also be derived from the variation of the energy functional 



Too r dr 



where Tqo is the local energy density 

1 



AttT, 



00 



(1-^2)2 

2r^ 



(7) 
(8) 

(9) 
(10) 
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Let me make two remarks which wih be useful in the subsequent discussion. First, note that, in 
general, the energy functional E is extremized only against variations with 64>{oo) = 0. However, 
it is also useful to consider more general variations for which 6(l){oo) ^ 0. Then, the variation 
of energy around a solution has the form 

SE = D (50(oo) , (11) 

where D = lim^^oo ^'^'P' is the dilaton charge. To avoid confusion I want to emphasize that the 
"surface term" on the right side of eq.(ll) is not of the Regge-Teitelboim type (in particular it 
cannot be cancelled by adding correction to energy) but it is rather a term which appears in 
variational problems with a free end. 

Second, note that the equations (7) and (8) have a "scaling" symmetry. Namely, if ■w{r) and 
0(r) are solutions so are 

w\{r) = w{e^r) , 

(12) 

(t>x{r) = 0(e^r) + A . 
Under this transformation the energy scales as follows 

E[wx,fx]=e-^E[w,f]. (13) 

The existence of this "scaling" symmetry excludes, via Derrick's argument, nontrivial static 
finite energy solutions with vanishing dilaton charge D (nota bene this also follows immediately 
from eq.(8)). However, when D ^ 0, Derrick's argument doesn't apply because for the variation 
induced by the transformation (12) (5^(oo) is nonzero, and therefore, as follows from (11), the 
energy is not extremized. Hereafter, I will assume that all solutions satisfy (f){oo) = 0, which 
can always be set by the transformation (12). This choice sets the scale of energy in the theory. 



3 Abelian solutions 

The equations (7) and (8) have two explicit abelian solutions. The first one is the vacuum 
solution 

w = ±l , (f) = (14) 

for which the energy has the global minimum E = 0. 
The second solution is 

w = , = ln(l + -) (15) 

r 

and its YM curvature is 

F = -Ts di9 AsinMif , (16) 

which corresponds to the Dirac magnetic monopolc with unit magnetic charge. There is also an 
electrically charged abelian solution related to (15) by the duality rotation : 

F = e~'^^ *F = T-^-^rg drAdt, 4> = -<l) ■ (17) 

(1 + ry 

These solutions have very interesting properties. The dilaton dramatically changes the properties 
of ^7(1) point monopoles (electric and magnetic). Without a dilaton, the energy density of a 
point monopole diverges at r = as Too ~ l/f^, whereas in the present case Too ~ 1/r^. Thus, 
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although the solution (15) is singular at r = 0, its total energy is finite and equals one! This 
result may be viewed as non-perturbative cancellation of two infinite self-energies: positive one 
of the point magnetic monopole and negative one of the dilaton. 

Since this solution will play an important role in the discussion of nonabelian solutions, it 
is useful to see how one can obtain it in a systematic way. Namely, in the w = sector, the 
energy (9) can be written as 

fOO 1 

E[(P] = / (r<^' + -e-<^)2 dr + e"^ \^ . (18) 
Jo r 

Thus, if (/)(0) = oo (and 0(oo) = 0), the energy is bounded from below by the value of magnetic 
charge (here set equal to one) and attains a global minimum = 1 on solutions of the first 
order Bogomol'nyi-type equation 

rV + e""^ = . (19) 

Solutions of this equation automatically satisfy the eq.(8) with w = 0. Elementary integration 
of this equation gives the solution (15). A detailed discussion of Bogomol'nyi inequalities in the 
Maxwell-dilaton theory, without an assumtion of spherical symmetry, will be given elsewhere 
[13]. 

4 Boundary conditions 

Now I will specify the boundary conditions for the globally regular solutions. They are deter- 
mined by the requirement that 

i) the local energy density be finite for all r 

Too < const < oo (20) 

which imposes boundary conditions for w and ^ at r = 0, and 

ii) the energy be finite 

E <oo (21) 

which imposes boundary conditions for w and (f> at infinity. 

It is easy to construct asymptotic solutions to eqs. (7) and (8) satisfying these boundary 
conditions. The solution near r = is 

w = l-br^ + 0{r^) , 

(22) 

= c - 26^2 + O(r^) . 
At r = oo the asymptotic solution is 

±w = 1 - d/r + 0{l/r^) , 

(23) 

(f) = e/r + 0{l/r^) . 

Here b, c, d, and e arc arbitrary constants. All higher order terms in the above expansions are 
uniquely determined, through recurence relations, by b and c in (22), and d and e in (23). 

Using these boundary conditions one can get the following elementary a priori results for 
the solutions of eqs. (7) and (8): 
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Lemma 1. The function w oscillates around zero between —1 and 1 (or \w\ = 1). 
Proof. It follows from eq.(7) that if w'{ro) = then at ro 

sgn w" = sgn w{w'^ — 1) . (24) 

This implies that w cannot have local maxima for w > 1 and local minima for w < —1. Since 
w{0) = 1 and 1^1^(00)1 = 1, this gives \w\ < 1 for all r > 0. Thus from (24), w"w < 0, which 
concludes the proof. 

Lemma 2. The function is monotonically decreasing. 
Proof. As above, this follows immediately from the maximum principle applied to eq.(8). 

Finally, note that, for the asymptotic behaviour (23), the radial magnetic curvature, B = 
Ts{l—w^)/r^, falls-off as and therefore all globally regular solutions have zero YM magnetic 
charge. 

5 Nonabelian solutions 

Let us assume that there exist a 2-parameter family of local solutions defined by the expansion 

(22) . Note that this is a nontrivial statement because the point r = is a singular point of 
the equations (7) and (8), so the formal power-series expansion (22) may have, in principle, a 
zero radius of convergence. A generic solution with initial data (22) certainly will not satisfy 
the asymptotic conditions (23) (in fact, the solution may even become singular at some finite 
distance). The standard numerical startegy, called the shooting method, is to find initial data 
(5, c) for which the local solution extends to a global solution with the asymptotic behaviour 

(23) . Actually, only 6 is a nontrivial shooting parameter, since one take arbitrary c and after 
finding the solution adjust the value of ^(oo) to zero using the transformation (12). For generic 
orbits with b < boo — 0.3795 the function w oscillates finite number of times in the region 
between w = —1 and w = 1 and then goes to ±00. For b > b^o all orbits become singular at a 
finite distance (in a sense that w' becomes infinite). 

The numerical results strongly indicate that there exist a countable family of initial data 
(bn:Cn), Ti ^ N, determining globally regular solutions Xn = {wn,4'n)- Here the index n labels 
the number of nodes of the function w. The values of initial data and energies of the first five 
solutions are given in Table 1. The functions w and (p for the first three solutions are graphed 
in Fig.l and 2. 

Table 1: Initial data (6, c) and energies of the first five globally regular solutions. 



n 


b 


c 


E 


1 


0.26083011 


1.711 


0.804 


2 


0.35351804 


3.374 


0.9659 


3 


0.3750017 


5.158 


0.9944 


4 


0.378754 


6.966 


0.9992 


5 


0.379373 


8.754 


0.99993 
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The solutions display three characteristic regions. The energy density Too is concentrated 
in the inner core region r < where i?i is approximately the location of the first zero of w. 
This region decreases with n and shrinks to zero as n — ^ oo. In the second region, Ri < r < R2, 
where R2 is approximately the location of the last but one zero of w, the function w slowly 
oscillates around w = with a very small amplitude. In this region the solution is very well 
approximated by the abelian magnetic monopole (15). This region extends to infinity as n — 00. 
Finally, in the asymptotic region r > R2, the function w goes monotonically to w = ±1 (hence 
the YM magnetic charge is gradually screened) and for r — 00 the solution tends to the vacuum 
{w = ±1,(P = 0). 

Because of these properties, the solutions are in striking resemblance to the BM solutions 
[6] of the Einsten-YM equations - the dilaton coupling has almost the same effect as the grav- 
itational coupling. In both cases the equilibrium configurations result from a balance between 
repulsive YM force and attractive, gravitational or dilatonic, force. There are indications that 
this analogy is even deeper. Below I will discuss two facets of the apparent duality between 
gravity and dilaton interacting with the YM field. 

First, I will show that in the YM-dilaton theory the energy of a static solution can be 
expressed as a surface integral at spatial infinity. To show this, I will first derive a simple virial 
identity. Consider a one-parameter family of field configurations defined by 



w/sir) = w{f3r) , 
ct,0{r) = <j>{l3r) . 
For this family the energy (9) is given by 

where 



1 f'^ 
2. Jo 



f 

Jo 



-24> 



dr . 



Since the energy is extremized at /3 = 1, it follows from (26) that on-shell 



Next, integrating eq.(8) one gets D = 
therefore eq.(29) yields 



(25) 

(26) 

(27) 
(28) 

(29) 



-2I2 {D is the dilaton charge defined in Section 2), and 
E = -D . (30) 



Thus the energy of a static solution can be read-off from the monopole term of the asymptotic 
expansion (23) of the dilaton field. This is a remarkable property which reminds very much the 
situation in general relativity and shows a relation between the dilaton field and the conformal 
degree of freedom of the metric. 

Secondly, the most striking analogy between our solutions and the BM solutions is their 
spectrum of energies (see Table 1 and compare with Table I in ref.[14]). In both cases the 
energies increase with n and are bounded from above by = 1. This cannot be a coincidence, 
but what distinguishes this particular value of energy which provides the common upper bound? 
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The answer is remarkably simple. The limiting Xqo solutions (whose energies give upper bounds) 
of our family and the BM family saturate the Bogomol'nyi inequalities in the abelian sectors of 
respective theories and therefore their energies are equal to the unit magnetic charge. To see this, 
consider first the dilatonic solutions. As was discussed above, the second region i?i < r < R2, 
covers the whole space as ra — > 00, since in this limit Ri ^ and R2 00. As n grows the 
amplitude of oscillations of the function w decreases and goes to zero as n — 00. Thus, for 
n — 00 the solution tends (nonuniformly) to the (singular) abelian magnetic monopole 
described in Section 3: 

X^ = {w = 0, (j) = ln(l + -) ) . (31) 

r 

As I have shown in Section 3, in the U{1) sector of the YM-dilaton theory, the static solutions 
satisfy the Bogomol'nyi inequality E > Q. The limiting solution X^^q saturates the bound in 
the Q=l sector. 

The behaviour of the BM solutions (in isotropic coordinates) is similar: as n — 00 the YM 
field tends to the abelian magnetic monopole while the metric develops a horizon and becomes 
the extremal Reissner-Nordstrom black hole solution with unit magnetic charge. Thus 

Xg^ = (w = 0, ds^ = -e-^^dt^ + e2^(dx2 + dy^ + dz'')) , (32) 

where 

U = ln(l + -) . (33) 
r 

It is well known that this solution saturates the Bogomol'nyi inequality in Einstein-Maxwell 
theory [15]. Actually, the limiting solutions (31) and (32) can be mapped one into another by 
the duality transformation U (f) and a 1/a in the abelian sector of the theory defined by 
the action (1). 

Prom the content of the last two paragraphs it is clear that to prove rigorously that the 
energy is bounded from above by one, one needs in the YM-dilaton theory (and in the Einstein- 
Yang-Mills theory in the case of BM solutions) a sort of Bogomol'nyi inequality with reversed 
sign, E < Q, which is saturated by the limiting abelian solution. Unfortunately, I wasn't able 
to find such an inequality. It would be probably easier, but also much less interesting, to find 
an upper bound which is not sharp (for the BM solutions that was done in ref.[14]). Also, it is 
not difficult to obtain not strict bounds on initial parameters. For example, multiplying eq.(8) 
by (f), integrating by parts and combining the result with eq.(29) yields the identity 



/•oo 

Jo 



2r2 



dr = 0, (34) 



which implies that (p{0) > 1. 



6 Stability analysis 

In this Section I address the issue of linear stability of the static solutions described above. To 
that purpose one has to analyse the time evolution of linear perturbations about the equilibrium 
configuration. I will assume that the time-dependent solutions remain spherically symmetric 
and the YM field stays within the ansatz (6). This is sufficient to demonstrate instability 
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because unstable modes appear already in this class of perturbations. The spherically symmetric 
evolution equations are 



1 



(35) 



,2^(1- W 



2\2 



0, 



(36) 



where dot denotes differentiation with respect to time t. 

Now, I take the perturbed fields: w{r) +dw{r,t), and ^(r) + (5^(r, t), where {w{r), f{r)) is a 
static solution, and insert them into the cqs. (35), (36). Linearizing and assuming the harmonic 
time-dependence for the perturbations: 6w{r,t) = e*'^*^(r) and 6(f){r,t) = e^'^*ip{r), one obtains 
an eigenvalue problem 



- e + 20'^' + 2w'i;' - ^(1 - 3w^)^ = (T^C , 



(r^V ) - 4e- 



-2(f) 



1mi-.^)C-(.'V^ 



2\2 



-)V 



2 2 / 



(37) 
(38) 



It is easy to check that if the perturbations satisfy the boundary conditions 



C(0) = 

e(cx)) =0 



^(0) = const , 

V^(oo) = , 



(39) 



then the above system is sclf-adjoint, hence eigenvalues are real. Instability manifests itself 
in the presence of at least one negative eigenvalue. 

To solve the eigenvalue equations (37), (38) with the boundary conditions (39), is a straighfor- 
ward but tedious numerical problem. I have done that for several lowest-energy static solutions. 
It turns out that the solution Xi has exactly one unstable mode of frequency o"^ ~ —0.0225. 
Each succesive static solution picks up one additional unstable mode (I have checked this up 
to n = 4). This is consistent with the fact that the limiting solution Xoo, given by (31), has 
infinitely many unstable modes. To prove this, consider the perturbations of with 5(f) = 0. 
Then, eq.(37) reads 



(40) 



r{l + r) 

This equation has infinitely many negative modes because the zero energy solution satisfying 
^(0) = 0, has infinitely many nodes as can be seen easily from the asymptotic solution. 

The result that the solution X„ has n unstable modes fits very nicely to the interpretation 
of solutions. In particular, for the interpretation of the solution Xi as a sphaleron, it is crucial 
that it has exactly one unstable mode. However, remember that I have considered a restricted 
class of perturbations, and by doing so some directions of instability might have been suprcssed. 
If there are additional unstable modes outside the ansatz (which I doubt), the interpretation of 
solutions given by Sudarsky and Wald would have to be revised. 
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7 Sudarsky and Wald's argument 

Sudarsky and Wald have recently proposed a heuristic argument which "explains" the existence 
and instability of the BM solutions of the Einstein- YM equations. This argument is, in principle, 
applicable to other theories involving the SU{2)-YM. field, which are not scale invariant and 
possess a stable solution. Below I outline the SW argument in application to the SU{2)-YM- 
dilaton theory. 

Let r be a space of all functions {Ai., <])), defined over R"' , for which the energy E is finite. Let 
r be a subspace of F with (/)(oo) = 0. The static solutions are extrcma of energy on T. One such 
extremum is the vacuum solution = 0) for which energy has a global minimum E = Q. 

Now, the key feature of the SU {2)-YM. group is the presence of "large gauge transformations" 
i.e. topologically inequivalent cross-sections of the YM-bundle, classified by the homotopy group 
T^siSU (2)) ~ Z. Thus, the energy functional E has a countable set of disconnected global minima 
corresponding to the trivial vacuum {Ai,(f) = 0) and all large gauge transformations of it. To 
avoid complications with the group of small gauge transformations G, it is convenient to pass 
from r to the space of gauge orbits F = T/G. 

Now, to apply Morse theory methods in Banach spaces, one needs a sort of compactness 
condition (like the Palais-Smale condition). A convenient way of implementing such a condition 
(which is here simply assumed to hold), is to introduce on F a Riemannian metric Gab (upper 
case latin letters denote indices of tensor fields on F), such that the flow generated by the vector 
field = -G^^VbE carries each point of F to a critical point of E. As discussed above, 
there exist a countable set of global minima of E. Since this set is disconnected, the flow defined 
by M'^ cannot carry all points of F to global minima (or local minima if any exist), because 
this would contradict the connectedness of F. Thus, the set, Fi C F, of points which do not 
fiow to local minima, must contain at least one critical point of E. A critical point on Fi with 
smallest energy Ei is a saddle point on F with exactly one unstable direction. This is believed 
to account for the existence of the solution Xi. 

Actually, there is a countable set of local minima of E restricted to Fi, namely Xi and all 
large gauge transformations of it. Hence, one can repeat the above argument, replacing F by Fi 
(and assuming that Fi is connected), to predict the existence of a submanifold F2 C Fi with 
a point X2, whose energy E2 minimizes E restricted to F2. The point X2 is an extremum of 
E on F, which has the 2-dimensional space of unstable directions. This is believed to account 
for the existence of the second static solution X2. All higher n solutions are predicted by the 
repetition of this argument. 

It seems very unlikely that the SW argument in its present form can be converted into a 
rigorous proof. However, the same argument can be made for spherically symmetric connections. 
Then, the powerful methods of equivariant Morse theory are available, and in fact these methods 
were succesfuUy applied in related problems [16]. In my opinion this is a very promising direction 
for future research. 

Another possibility of proving rigorously the existence of numerical solutions found in this 
paper is to apply the methods of the dynamical systems theory. This approach was used recently 
by Smoller and his collaborators [17,18], who proved the existence of the BM family of solutions 
to the Einstein- YM equations. A similar proof should be possible for the YM-dilaton equations 
although it might be more difficult, because here the corresponding (nonautonomous) dynamical 
system is four-dimensional whereas in the Einstein- YM case it is three-dimensional. 
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Figure captions 

Fig.l The function w for the first three globally regular solutions. 
Fig. 2 The function (p for the first three globally regular solutions. 



